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Abstract 

The modern nucleon-nucleon interaction of Myers and Swiatecki, 
adjusted to the properties of finite nuclei, the parameters of the mass 
formula, and the behavior of the optical potential is used to calculate 
the properties of /3-equilibrated neutron star matter, and to study the 
impact of this equation of state on the properties of (rapidly rotating) 
neutron stars and their cooling behavior. The results are in excellent 
agreement with the outcome of calculations performed for a broad 
collection of sophisticated nonrelativistic as well as relativistic models 
for the equation of state. 
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I Introduction 



The equation of state (EOS) of dense nuclear matter plays a decisive role 
in many branches of physics. It is an important input quantity for the un- 
derstanding of heavy-ion collisions, supernova explosions, and the structure 
of neutron stars [H-lH. Therefore great efforts have been made by several 
groups to derive the EOS of such matter. These models can be divided into 
two categories, namely nonrelativistic potential models and relativistic, field- 
theoretical ones. A frequently used model that belongs to the first category is 
the Skyrme Hamiltonian which has the advantage, being a phenomenological 
approach, to reproduce the nuclear data rather well. Extensive calculations 
of neutron star matter properties for this scheme have been performed (see, 
for instance, Refs. [^, H, |^). More elaborate, microscopic schemes constitute 
the nonrelativistic potential models, where one uses realistic nucleon-nucleon 
potentials in combination with sophisticated many-body approximations. A 
typical example of this kind is the work of Wiringa et al. p], in which the 
properties of matter were calculated for a Hamiltonian containing two- and 
three-body interactions. These extensive calculations with some later refine- 
ments are still the basis for the calculation of neutron star properties (in the 
nonrelativistic approach) in several investigations |0. The disadvantage of 
such a microscopic treatment is the numerical complexity of the method. For 
that reason it is very tempting to use simpler models which are easier to deal 
with, and make comparisons with respect to the properties of finite nuclei, 
the parameters of the mass formula, and neutron stars. For this purpose we 
selected the new Thomas-Fermi approach (TF) of Myers and Swiatecki |jTO| , 
where the density- and momentum dependent interaction is given by [the 
upper (lower) sign corresponds to nucleons with equal (unequal) isospin]: 

V..r = -^/f^) (1.1) 




Po V a 

The quantities Po, Po, and Tq^{= pl/2m^) denote the baryon number density, 
Fermi momentum, and the kinetic single-particle energy of symmetric nu- 
clear matter at saturation, respectively. The potential's radial dependence, 
/, is chosen to be Yukawa type. The choice ^ 7^ C leads to a better description 
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of asymmetric nuclear systems, and the behavior of the optical potentials is 
improved by the term a{2p/ po)"^/^ , where p^/^ = 0.5(pi^^+p2^^)- As discussed 
in great detail by Myers and Swiatecki, this model (with only seven free pa- 
rameters) has the advantage over the former Seyler-Blanchard interaction 
|TT| , to not only reproduce the ground-state properties of finite nuclei 
and symmetric infinite nuclear matter but also the optical potential and, as 
revealed by a comparison with the theoretical investigations of Friedman and 
Pandharipande [|l3l, the properties of pure neutron matter, too. Therefore 
we think this approximation, which is able to account for so many features in 
a very satisfying fashion, constitutes, despite its simplicity, an excellent can- 
didate for the investigation of dense matter in the nonrelativistic approach 
(see also where a simpler force has been used). Since the standard com- 
parison with nuclear physics quantities has already been performed by Myers 



and Swiatecki |jTO[ , we will concentrate in this contribution on the properties 
of neutron star matter and the structure and thermal evolution of neutron 
stars constructed for their model. The outcome will be compared with the 
one obtained for a broad collection of other models for the equation of state, 
which comprises nonrelativistic potential models as well as relativistic, field- 
theoretical ones. We shall not go into details concerning these models, whose 
properties have been discussed in numerous investigations prior to this one 
(see, for instance, Refs. §, |, i and [|15l-[gTl). 



II Equation of State 

a Symmetric and Asymmetric Nuclear Matter 

The explicit expressions for the energy per baryon, E/A, chemical potentials, 
p, pressure, P, effective masses, m*, and symmetry energy J are given in 
Refs. ||I0| , p2| and will not be repeated here again. For the parameters of the 
TF model of Myers and Swiatecki we use the most recent set denoted TF 
p3| (We denote the different parameter sets with TF90, TF94, and TF96. 
If not stated otherwise the latest parametrization, TF96 (= TF), is used.). 
Earlier parametrizations (TF90, TF94) have larger statistical errors in fitting 



the nuclear data. The EOSs for the parameter sets of Ref. are rather 
similar, but the EOS for the parameter set TF90 of Ref. is stiffer than 
the latter EOSs (see Table I). For the purpose of comparison, we give the 
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parameters of symmetric infinite nuclear matter for the different models in 
Table I. (An overview of all EOSs used in this work is given in Table II). 

Nonrelativistic Brueckner-Hartree-Fock calculations (BHF) with two- 
body forces are in general not capable to reproduce the correct saturation 
properties. The resulting EOSs are too soft. This deficiency is corrected by 
adding three-body forces (TBF) which lead to improved saturation proper- 
ties and a stiffer EOS. Nevertheless the BHF calculations occur to be not 
satisfactory enough for the description of symmetric infinite nuclear matter 
(see Refs. p|, |^ and Table I) and we will therefore not discuss them in greater 
detail. Stiffer EOSs are also obtained in the Dirac-Brueckner-Hartree-Fock 
approach (DBHF) ||2l[],p4|-[|28l. The general trend seems to be that varia- 
tional calculations, BHF supplemented with TBF, and DBHF calculations 
agree with each other more or less only up to moderate densities. At higher 
densities one obtains a closer agreement of DBHF with the variational cal- 



culations, depending on the matter's asymmetry [g, ^ ^ Our main 
intention is to discuss the properties of the EOS of neutron star matter, and 
compare the results with those obtained in the nonrelativistic many-body 
approach as well as the fieldtheoretical treatment. As examples belonging to 
the first category, we have selected EOSs from the investigations of Wiringa 
et al. 0, which have been used frequently in neutron star calculations by 
others, too. The two samples are (cf. Table II) the variationally-based EOSs 
computed for the Urbana UV14 two-nucleon potential supplemented with 
the UVII three-nucleon potential (WUU), and the density-dependent UV14 
potential supplemented with the three-nucleon interaction TNI (WUT). As 
representatives of the DBHF and relativistic Hartree-Fock (RHF) scheme we 
have chosen EOSs from Refs. [T|, |2l|, ||, |6|. 

As pointed out above, the DBHF-EOSs are stiffer than BHF-EOSs based 
on two-body forces, but by introducing three-body forces one can simulate 
the relativistic behavior to a certain extent p|, § . In Fig. 1 we illustrate the 
situation for symmetric nuclear matter. The agreement of all models, with 
the exception of the BHF model based on two-body forces 0, is 
better than in the case of pure neutron matter (see Figs. 2 and 3), where 
larger deviations occur above p = 0.4 fm~^ . The potential model of Myers 
and Swiatecki leads to a pure neutron-matter EOS which is located between 
the two nonrelativistic variational calculations WUU and WUT. The larger 
differences between nonrelativistic and relativistic models for higher asym- 
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metry are explainable by the different behavior of the symmetry energy. The 
nonrelativistic potential models cause a stronger repulsion in isospin singlet 
states for higher densities which makes for increasing asymmetry the 
nonrelativistic EOSs softer in comparison with the relativistic ones. This 
behavior can be clearly seen in the density dependence of the symmetry en- 
ergy (see Fig. 4), where TF shows a similar behavior as WUT. It seems to 
be a general feature that relativistic treatments lead to a monotonic increase 
of the symmetry energy, in contrast to the nonrelativistic models where the 
symmetry energy saturates (for BHF H, ||, |^), or even decreases with in- 
creasing density (variational calculations). As we shall discuss later, this 
feature has a severe impact on the composition of neutron star matter, since 
the symmetry energy determines the proton fraction of neutron star matter. 
The proton fraction in turn plays an important role for the cooling behavior 
of neutron stars, since proton fractions above ~ 0.11 - 0.13 permit stars to 



cool very efficiently via the direct Urea process |29, 30 



b Neutron Star Matter 

Since neutron stars are bound by gravity, which is much weaker than the 
Coulomb force, and have life times practically infinite compared to the char- 
acteristic weak interaction time scale, ~ 10~^° s, neutron star matter is 
subject to the constraints of charge neutrality and generalized /3-equilibrium, 
respectively |, g g 0, 



^ ,,(2J, + 1)^- 2 |i = 0. (11.2) 

B=p,n L=e,fi 

and 

jJ'B = jJ'n- QB /iL , ^lf,= He ■ (11.3) 

Here Hb^Ib) denotes the chemical potentials (electric charges) of baryons, 
and /iL the chemical potentials of leptons e~, Due to these constraints 
the calculation of neutron-star-matter properties differs from the treatment 
of asymmetric matter, since the composition has to be determined selfcon- 
sistently subject to the two additional constraints ( |1I.2|) and ( |1I.3|) . The new 
baryon/lepton degrees of freedom lower the energy and pressure of neutron 



star matter in comparison with pure neutron matter y, |18|, |2^, |3T 
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This feature is illustrated in Fig. 5, where the pressure of neutron star mat- 
ter is shown as function of density and compared with the pressure of pure 
neutron matter. Due to the new degrees of freedom, the pressure is lower in 
/3-stable matter. To demonstrate the influence of muons on the EOS, we also 
show the EOS computed for n, p, and e~ only. It indicates that the muons 
play only a minor role for the stiffness of the EOS. Nevertheless they are 
important for the composition. For the neutron star matter EOS computed 
for the model of Myers and Swiatecki one does expect deviations from the 
microscopic calculations which are of the same order as for the calculations 
performed for fixed asymmetry. 

In Figs. 6 and 7 we show the energy per baryon of neutron star matter for 
the different models. The nonrelativistic models give in general softer EOSs, 
but the model of Myers and Swiatecki comes rather close to the relativistic 
results (HWWl), especially at lower densities. The earlier parametrization 
of Myers and Swiatecki, that is TF 90, comes even closer to the relativistic 
outcome since the EOS is much stiffer. 

The same behavior holds also for the pressure, which is exhibited in Fig. 8. 
Much larger differences than for the EOS occur for the composition of neu- 



tron stars 1^ . Due to the behavior of the symmetry energy discussed above, 
one encounters in relativistic models a monotonous increase of the proton 
fraction, x (= Pp/p) with density (in the case of no hyperons and/or meson 
condensates), in contrast to nonrelativistic microscopic theories and the TF 
model of Myers and Swiatecki for which x decreases (or saturates as for BHF) 
at higher densities. This is illustrated in Fig. 9 and Table III, where the com- 
position is compared for several different models. In this context we note that 
the stiffness of modern fieldtheoretical models which - besides neutrons and 
protons account for the population of more massive baryon states (especially 
hyperons) too - comes closer to the stiffness one of nonrelativistic neutron 
star matter EOSs [0, |19|, 0. Here one obtains again a decrease of the 
electron/muon fraction, since it is energetically favorable for the system to 
achieve charge neutrality among the hyperons themselves. 



Ill Neutron Star Properties 

In order to determine neutron star properties one has to solve Einstein's 
equations, for which knowledge of the EOS, i.e. -P(e), is necessary. Einstein's 
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curvature tensor G^i, {R^u, Ofiu, R, and T^,^ denote the Ricci tensor, metric 
tensor, Ricci scalar, and energy-momentum tensor density, respectively) is 
determined by 

G^. = R^.u-\ 9f.u R = SvrT^, (e, P(e)) . (III.4) 

For the determination of the properties of (rapidly) rotating neutron stars 
and their limiting periods, one has to generalize the Schwarzschild metric of 
a spherically symmetric, static star to the one of a rotationally deformed, 
axially symmetric body |53[ : 



Q denotes the neutron star's rotational frequency, and ci;(r, 6; Q) is the angular 
velocity of the local inertial frames (dragging of the local inertial frames) 
1,11,11). 

Of great interest with respect to the identification of rapidly spinning pul- 
rotating neutron stars is the maximum possible rotational frequency 
of such an object, for which no simple stability criteria exist in general rel- 
ativity. However an absolute upper limit is given by the Kepler frequency, 
Qk, at which mass sedding at the equator sets in. It is determined as the 
solution of the following transcendental equation [^, [S^, Q , 



L J eq 

which is to be evaluated at the star's equator. V denotes the velocity of 
a particle rotating at the star's equator. Equations (|III.4|) and ( [1II.6| ) are 
to be solved simultaneously by a selfconsistent numerical iteration scheme, 
since neither the metric functions in ( |III.5| ) nor the frequencies fix and u (all 
depend on the star's unknown mass) are known. 

In Fig. 10 we exhibit the neutron star masses as a function of central 
energy density for different EOSs. Since WUU and WUT behave rather 
similar up to mass densities of about ~ lO^^g/cm^, the mass-density curves 
computed for these two EOSs are very similar to each another. The devi- 
ations beyond that density arise from the different three-body forces. The 
TF-EOS behaves for smaller energy densities like the relativistic EOS, but 
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shows at higher densities features similar to a medium-stiff EOS. So the 
properties of neutron stars computed for TF and HWW2, which are hsted 
in Tables IV and V, are similar to each other as long as the star's mass is 
below about 1.4 Mq. Besides HWW2 the EOS of Brockmann and Machleidt 
P^ , |2B| has been chosen here as a further representative of a relativistic EOS. 
The former accounts for the presence of hyperons in neutron star matter too 



27| , |31[| . The stiffness of an EOS at high densities is known to play a key role 
for the maximum possible mass that a neutron star can have. Since the non- 
relativistic EOSs tend to be less stiff at high densities than their relativistic 
counterparts (the only exception is WUU), we obtain less massive neutron 
stars for the nonrelativistic EOSs. 

Another effect caused by the stiffness of an EOS concerns the radius 
of the stellar configurations. In general, the stiffer the EOS the lower the 
star's central mass density and thus the lower the gravitational force that 
pulls the star together. So stars constructed for stiffer EOSs will have larger 
radii than those constructed for softer EOSs. This behavior is confirmed 
in Fig. 11, where the radius-mass relations of neutron stars constructed for 
TF and HWW2 are shown: the stiffer model for the EOS, that is HWW2, 
leads to somewhat bigger stars, provided their mass is beyond about one 
solar mass. The masses of stars lighter than this value are determined by 
the stiffness/ softness of the EOS at nuclear and intermediate densities where 
TF behaves stiffer than HWW2, leading to neutron stars that are somewhat 
bigger for TF than for HWW2. The radial changes due to rotation at Qk 
amount about 2 km, except for the very light stars of each sequence. The 
mass increase of the heaviest stars due to rapid rotation is about 0.3 Mq (see 
also Fig. 12). Various properties of static as well as rotating neutron stars, 
computed for a number of different EOSs, are listed in Tables VI- VIII. All 
properties refer to neutron stars of mass M = 1.4 Mq (static and rotating), 
about which the observed masses tend to scatter. 

In Fig. 13 we show the density profiles of nonrotating neutron stars of 
mass M = 1.4 Mq, which reflect what has just been discussed above in 
connection with the dependence of stellar radii on the stiffness/softness of 
the EOS: The microscopic EOSs WUU and WUT, behaving somewhat softer 
than TF and HWW2, lead to smaller stars with a higher central density (for 
more phenomenological relativistic EOSs, see Ref. [^). Moreover, as can 



be inferred from Fig. 13, the TF-EOS density profile is rather close to the 
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one of the relativistic HWW2-E0S because of the rather similar behavior of 
these EOSs over the energy densities relevant for a star of this mass. 

A final remark concerns the velocity of sound, which does not exeed the 
causal limit in the TF model for the range of relevant neutron star densities, 
i.e., Cs < 0.96 c for p < 1.246 fm~^. The causal limit is reached at p = 
1.36 fm~^, which is considerably higher than the highest density reached in 
the most massive star of the sequence. 

As mentioned above, the Kepler frequency gives an absolute upper limit 
on the critical rotational frequency of a neutron star. Its rotation may also 
be limited by the gravitational radiation-reaction instability. Since the the- 
ory of this effect is too lengthy, we will not outline the treatment in this 
investigation but refer to 0, ^ |36| for details. In Fig. 14 we exhibit the lim- 
iting rotational periods set by both the Kepler criterion and the gravitational 
radiation-reaction instability. The latter depend on temperature, for which 
we have chosen representative values of T = 10^ K, corresponding to an old 
neutron star, and 10^*^ K typical for a young newly formed one. The limiting 
periods set by the gravitation-radiation instability are larger than the Ke- 
pler periods and thus set a more stringent limit on rapid rotation than mass 
shedding. Moreover they increase with temperature, since viscosity damps 
the instability modes less efficiently in hot stars. A comparison between TF 
and HWW2 shows that the differences in rotational periods are rather small, 
except for hot neutron stars where somewhat larger deviations show up. The 
shaded rectangle in Fig. 14 covers masses and periods of observed pulsars 



||37|| . Hence both EOSs, based on the assumption that neutron star matter 
is made up of hadrons and leptons rather than other kinds of exotic forms of 
matter 0] , are in accordance with the body of presently existing data. 

A further important facet of neutron star physics concerns the cooling 
behavior of such objects. In Ref. we have already calculated and dis- 
cussed in great detail their cooling behavior for different EOSs depending 
on the involved processes (see, for instance. Table 5 of Ref. [^). Here we 



show and compare the cooling tracks of different neutron star models of mass 
M = 1.4 Mq. The underlying EOSs are the nonrelativistic TF and WUU 
models, and the relativistic HWW2. As an example of an enhanced neutrino 
emission process we consider the direct Urea process, which is only possi- 
ble if the proton fraction exceeds a certain critical value (~ 0.13, see Refs. 



29|, |30|). Among the considered EOSs, only HWW2 allows for the direct 
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Urea process. This leads to a large drop of the star's surface temperature 
at r ^ 30 yr (see the dashed curve of the nonsuperfiuid model in Fig. 15), 
since the core cools down very fast via the enhanced neutrino emission pro- 
cess in comparison with the crust, causing a temperature inversion in young 
stars. Depending on the crust thickness, the cooling wave formed by the tem- 
perature gradient reaches the surface and causes the sharp decrease of the 
surface temperature, followed by a rather flat behavior of the cooling tracks 
up to 10'' years ||35|, We show also the behavior of superfluid neutron 
stars. Superfluidity reduces the neutrino emissivity, the heat capacity, and 
the thermal conductivity by an exponential factor of exp(— A/fcT), where A 
denotes the gap energy (see Table 4 of Ref. ||3^ for the used gap energies). 
The observational data are described in Ref. (see Table IX). The 



obtained effective surface temperature depends crucially on whether a hy- 
drogen atmosphere is used or not. Since the photon flux, measured solely in 
the X-ray energy band, does not allow to determine what atmosphere one 
should use, we consider both the blackbody model (sohd error bars in Fig. 15) 
and the hydrogen-atmosphere model (dashed error bars). The plotted errors 
represent the Sa error range due to the small photon fluxes. The pulsars 
ages are determined by their spin-down times assuming a canonical value of 
3 for the breaking index. In reality the breaking index may be quite different 
from 3. Its variation between 2 and 4, for instance, would change the age 
of Geminga as indicated by the horizontal error bar shown at the bottom of 
Fig. 15. 

The cooling tracks of the nonrelativistic models, that is TF and WUU, are 
almost identical. Since the relativistic HWW2 model cools mainly through 
the direct Urea process in the neutrino cooling era, its cooling behavior differs 
considerably from the nonrelativistic models. Superfluidity reduces the neu- 
trino emissivity which leads to a higher surface temperature in the neutrino 
cooling era. Later, in the photon cooling era, the surface temperature be- 
comes smaller for the superpfluid models. The observational data can almost 
perfectly be described by both the nonsuperpfluid and the superfluid non- 
relativistic models as well as the superfluid relativistic model, provided one 
assumes that the pulsars have no hydrogen atmosphere (except PSR 1055- 
52, which could be explained by internal heating; see, e.g., Refs. p]] ). 
However, if some of the pulsars prove to have a hydrogen atmosphere, these 
models seem to be too hot. In this case superfluid enhanced cooling or inter- 
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mediate cooling might be necessary. Whether the observed pulsars have 
a hydrogen atmosphere could be decided if one considers multiwavelength 
observations, as suggested by Pavlov et al. P2 |. 



IV DISCUSSION AND CONCLUSIONS 

We have calculated within the modern nonrelativistic TF model of Myers 
and Swiatecki the properties of neutron star matter and compared the out- 
come with microscopic nonrelativistic and relativistic treatments. The main 
purpose of this contribution was to test this model with respect to its im- 
plications for the physics of neutron stars. It turned out that, despite its 
simplicity, the modern TF model shows the same features with respect to 
the EOS of neutron star matter as complicated nonrelativistic variational 
calculations. Moreover the model is in even closer, over the relevant density 
regions encountered in neutron stars, with microscopic relativistic EOSs. The 
composition of neutron star matter obtained for TF is similar to other nonrel- 
ativistic models for which the symmetry energy saturates or even decreases at 
higher densities too, favoring so, in contrast to relativistic models, again an 
increase of the neutron fraction at higher densities. The resulting gross struc- 
tural properties of neutron stars (masses, redshifts etc.) are in accordance 
with observations. Also interesting is the finding that the limiting rotational 
neutron star periods computed for TF, set by either the Kepler criterion or 
the gravitational-radiation reaction instability, agree rather closely with the 
results of modern relativistic EOSs. The limiting period of a M ~ 1.4 Mq 
neutron star is about 1 ms, which too is in agreement with the present ob- 
servations. Limiting periods of this size seem to be a general feature of EOSs 
computed hadronic/leptonic matter |P, |2^. Finally we studied the thermal 



evolution of neutron stars for TF. Here, somehow surprisingly, it turned out 
that no enhanced cooling mechanism needs to be invoked to achieve agree- 
ment between theory and observation; cooling via only the modified Urea 
process occurs to be sufficient. (We recall that cooling via the direct Urea 
process is not permitted for TF, as is the case for the nonrelativistic micro- 
scopic neutron star matter EOSs stdied in this work.) 

From these findings we conclude that the modern Thomas-Fermi model 
of Myers and Swiatecki is not only very suitable for the description of low- 
energy nuclear physics but for phenomena in high-density regime too! Our 
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study reveals that the TF-EOS of neutron star matter shows similar features 
as sophisticated microscopic nonrelativistic EOSs and, over some density 
range, comes even close to microscopic relativistic EOSs. The neutron star 
properties and the cooling behavior derived for TF are in excellent agreement 
with the body of presently existing data. Due to these surprising features we 
think that the model of Myers and Swiatecki is also a very good candidate for 
the nonrelativistic (many-body) treatment of neutron stars and, therefore, 
should be included in further nonrelativistic investigations in this field. The 
handling of this model with only seven parameters is relatively simple in 
comparison with elaborated microscopic schemes, which makes it an ideal 
candidate for this kind of astrophysical studies. 
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Table captions 



TABLE I: Saturation density, poo, energy per baryon, E /A, incompressibil- 
ity, K, and symmetry energy, J, of symmetric infinite nuclear matter 
for different models. These are TF: Thomas-Fermi model of Myers and 



Swiatecki fl^, |2^; WUU and WUT: variational calculation performed 
for UV14 plus UVII and UV14 plus TNI, respectively g |; BHF2 and 
BHF3: Brueckner-Bethe approximation with two- (AV14) and three- 
body forces (AV14+UVII), respectively % |]; DBHF^^) and DBHF^^); 
relativistic Dirac-Brueckner-Hartree-Fock approximation with solutions 



in the space of positive energy-spinors only and the full Dirac space 



21[] , respectively (A and B denote different Brockmann-Machleidt po- 
tentials pT| , p4[); NLl and NL-SH: phenomenological relativistic mean 
field theory for the parameter sets of Reinhard and Sharma, respec- 



tively 0); SKM* and SIII: Skyrme forces ||. 

TABLE II: Equations of state (EOSs) used in this work. 

TABLE III: Proton fraction, a;, and energy per baryon (in MeV) of beta- 
stable matter (neutrons, protons, electrons and muons). 

TABLE IV: Neutron star properties for beta-stable TF model (neutrons, 
protons, electrons and muons). Mci^ec) denotes the gravitational mass 
in solar mass units as a function of central mass density. The central 
star pressure is denoted pc- The amu mass Ma minus the gravitational 
mass Mg is effectively the binding energy liberated when the NS is 
formed. R stands for the star's radius, denotes the stellar crust using 
2.4x10^^ gcm~^ as the boundary, / denotes the moment of inertia, 
and z the surface redshift. For a comparison with WUU and WUT see 
Ref. i. 

TABLE V: Same as Table IV, but for the HWW2 model. 

TABLE VI: Comparison of the properties of a static, spherical NS of mass 
1.4 Mq for different models. 

TABLE VII: Properties of rotating neutron star models of mass M = 
1.4 Mq and rotational period P = 1 ms, calculated for different EOSs. 
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The entries are (from top to bottom): central energy density, ec'-, per- 
cental mass increase relative to nonrotating star model of same 
AM/M\ equatorial and polar radii, i?eq and i?p, respectively; moment 
of inertia, /; stability parameter, t; injection energy, /3; redshift at 
the pole, eccentricity, e; quadrupole moment, 11. The influence of 
different parametrizations of the TF model on NS properties is shown 
too. The stiffness of the TF EOSs decreases from TF90n (pure neutron 
matter, represents the stiffest EOS) to TF96. 

TABLE VIII: Same as Table VII, but for a rotational period of P = 1.6 ms. 

TABLE IX: Surface temperatures, Tj^, and spin-down ages, t, of several 



observed pulsars pq, 38, 39 . 
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TABLE I 



Method 




E/A 


Poo 


K 


J 






(MeV) 


(fm-3) 


(MeV) 


(MeV) 


TF90 




-16.53 


0.165 


301 


32 


TF94 




-16.04 


0.161 


234 


32 


TF96 




-16.24 


0.161 


234 


33 


WUU 




-15.5 


0.175 


202 


30 


WUT 




-16.6 


0.157 


261 


29 


BHF2 




-17.8 


0.280 


247 




BHF3 




-15.2 


0.194 


209 




DBHFA(i) 


-15.6 


0.185 


290 




DBHFA(2) 


-16.5 


0.174 


280 


34 


dbhfbw 


-13.6 


0.174 


249 




DBHFB(2) 


-15.7 


0.172 


249 


33 


NLl 




-16.4 


0.152 


212 


43 


NL-SH 




-16.3 


0.146 


356 


36 


SkM* 




-15.8 


0.160 


216 


30 


SIII 




-15.9 


0.145 


355 


28 


a) in Ref. 




-15.7 


0.158 


246 
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TABLE II 



EOS 


Till" PTi^ cf^c^Y\ 

±11 LCI CLl^ LlVjll 


ivj-cxiiy uwu.y 


Ref. 










TF90 


A/Tvprc; anr] Swii^tpplci 

J.VJ. y VjX o cajIJ.*^ kj vv icjl \j\jK^±\.i 


Tb nm P^ ^-Eprrn i 


|10 




TPQzL 


iV±y t;i b cLiiU. O WlcLLcL^ivl 


± iioiiiao--r Cl iiii 


|23 






Ix/T'xrPTc; Q n H ^^sn n "I'Pr'ln 

iVl y Cl O dllLl kJ W ICL tCL/lvl 


1. livjiiido X Cl liii 


P 




Willi 


T It iisi Ti V—Tii 1 r*! pr^Ti T^^~»"f"PTl"^ 1 si 1 T T \/ \ A. T\\^^Q. 

KJ L UcXiLLcXi llUL-lv/Wll UVjljClluldl \J V Xt: UlLlo 


vdi idtivJiicLi 


i 




\JLUoiLLcl O ilLIL-lcUii UULt;liLlcLl KJ Vll 






WITT 


UiUdiid Zi-iiLIL/lcUii UULt;liLldl KJ V Xt: 


VcLl IcLLlUlldl 


i 




Tlrl^ririri .'^-Tinrlpori nntpntiril TTsFT 






RHF9 


A TcroTi Ti P V—Tin r'lpoTi T»o"i'P'n"i"i fi 1 A \/ 1 A 

X1.1 tivjllllt:!' Z( llLll^lCvJll [J Vj LClluldl i V V Xt: 


llvjlll CldulV loljlL/ 


i 


i 






Brueckner-Bethe 




BHF3 


Argonne 2-nucleon potential AV14 plus 


nonrelativistic 


B 


i 




Urbana 3-nucleon potential UVII 


Brueckner-Bethe 




RE 


Brockmann-Machleidt potential A 


relativistic 


24 


M 






Brueckner- 








Hartree-Fock 




HWWl 


Brockmann-Machleidt potential A 


relativistic 


[28 






(no hyperons) 


Brueckner- 








Hartree-Fock 




HWW2 


Brockmann-Machleidt potential B 


relativistic 


27 






+ HFV 


Brueckner- 








Hartree-Fock 




HV 


Exchange of a, and 


relativistic 


i 




p mesons (with hyperons) 


Hartree 




HFV 


Exchange of a, u; and 


relativistic 


i 




p mesons (with hyperons) 


Hartree-Fock 
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TABLE III 



p (im j 


TF 


WUU 


WUT 


x{p) 


E[p,x) 


x[p) 


E{p, x) 


x{p) 


E{p, x) 


U.U7 


0.0276 


5.88 


0.0193 


8.13 


0.0247 


6.25 


0.08 


0.0308 


6.25 


0.0213 


8.66 


0.0273 


6.43 


0.10 


0.0365 


7.05 


0.0253 


9.70 


0.0317 


6.87 


0.125 


0.0427 


8.25 


0.0300 


11.06 


0.0357 


7.69 


0.15 


0.0502 


9.48 


0.0345 


12.59 


0.0388 


8.89 


0.175 


0.0572 


10.97 


0.0402 


14.18 


0.0418 


10.32 


0.20 


0.0630 


12.83 


0.0464 


15.92 


0.0442 


12.13 


0.25 


0.0715 


17.78 


0.0572 


20.25 


0.0469 


16.76 


0.30 


0.0767 


24.36 


0.0632 


25.78 


0.0476 


22.53 


0.35 


0.0794 


32.57 


0.0673 


32.60 


0.0476 


29.18 


0.40 


0.0803 


42.34 


0.0741 


40.72 


0.0472 


36.75 


0.50 


0.0780 


66.28 


0.0854 


61.95 


0.0401 


56.06 


0.60 


0.0725 


95.56 


0.0959 


90.20 


0.0311 


79.19 


0.70 


0.0651 


129.62 


0.1108 


126.20 


0.0195 


106.04 


0.80 


0.0572 


167.97 


0.1215 


170.50 


0.0054 


135.46 


1.00 


0.0417 


256.15 


0.1239 


291.10 


4.8e-4 


200.89 


1.25 


0.0274 


385.57 


0.0855 


501.00 


0.0000 


294.00 


1.50 


0.0152 


534.64 


0.0195 


753.00 


0.0000 


393.00 
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TABLE IV 



(1014) 
(g/cm^) 


Pc 

(10^4) 
(dyn/cm^) 


Mg 
(Me) 


Ma-Mg 
(Me) 


R 
(km) 


(km) 


/ 

(10^4) 

(g cm^) 


z 


2.5 


0.2703 


0.131 


0.001 


17.89 


16.08 


0.569 


0.011 


3.0 


0.4496 


0.184 


0.003 


14.17 


9.90 


0.854 


0.020 


3.5 


0.7009 


0.248 


0.005 


12.68 


7.06 


1.26 


0.030 


4.0 


1.046 


0.327 


0.009 


11.99 


5.36 


1.82 


0.043 


5.0 


2.007 


0.504 


0.022 


11.53 


3.54 


3.30 


0.072 


6.0 


3.387 


0.696 


0.041 


11.45 


2.57 


5.15 


0.104 


7.0 


5.218 


0.887 


0.067 


11.45 


1.98 


7.21 


0.139 


8.0 


7.518 


1.07 


0.098 


11.45 


1.59 


9.26 


0.174 


10.0 


13.46 


1.36 


0.163 


11.38 


1.13 


12.9 


0.244 


12.5 


23.28 


1.62 


0.239 


11.19 


0.82 


16.0 


0.322 


15.0 


35.38 


1.78 


0.297 


10.95 


0.64 


17.9 


0.387 


17.5 


49.47 


1.88 


0.339 


10.71 


0.53 


18.7 


0.442 


20.0 


65.27 


1.94 


0.366 


10.47 


0.46 


18.9 


0.487 


25.0 


101.1 


1.99 


0.392 


10.06 


0.37 


18.6 


0.553 


30.0 


141.6 


2.00 


0.396 


9.73 


0.33 


17.9 


0.598 


35.0 


185.6 


1.99 


0.390 


9.43 


0.27 


17.0 


0.633 
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TABLE V 



fg/cm'^) 

Vo/ ^^^^ 1 


Pc 

fdvn/cm^) 


Mr 
(Mn) 


Ma - Mr 
(Mn) 


R 
(km) 


(km) 


/ 

(10^4) 

vo ^^^^ ) 


z 


2.5 


0.3146 


0.144 


0.002 


14.47 


12.44 


0.598 


0.015 


3.0 


0.5156 


0.212 


0.004 


12.59 


8.03 


1.04 


0.026 


3.5 


0.7731 


0.282 


0.007 


11.95 


6.11 


1.54 


0.037 


4.0 


1.055 


0.345 


0.010 


11.68 


5.06 


2.02 


0.047 


5.0 


2.218 


0.539 


0.024 


11.37 


3.29 


3.66 


0.078 


6.0 


3.655 


0.735 


0.048 


11.36 


2.41 


5.59 


0.112 


7.0 


5.724 


0.955 


0.085 


11.44 


1.83 


8.09 


0.152 


8.0 


9.006 


1.21 


0.146 


11.52 


1.37 


11.4 


0.204 


10.0 


17.48 


1.61 


0.280 


11.55 


0.91 


17.0 


0.302 


12.5 


30.81 


1.91 


0.423 


11.42 


0.64 


21.5 


0.404 


15.0 


46.21 


2.06 


0.519 


11.20 


0.50 


23.6 


0.481 


17.5 


61.29 


2.13 


0.571 


11.00 


0.43 


24.2 


0.530 


20.0 


78.08 


2.17 


0.602 


10.81 


0.40 


24.2 


0.568 


25.0 


112.3 


2.19 


0.624 


10.45 


0.33 


23.4 


0.619 


30.0 


148.3 


2.18 


0.621 


10.18 


0.31 


22.2 


0.647 


35.0 


185.1 


2.15 


0.609 


9.95 


0.29 


21.1 


0.665 



TABLE VI 



Quantity 


TF 


wuu 


WUT 


HWW2 


6^(10^4 g/cni3) 


10.30 


10.42 


12.12 


8.9 


Pc(10^^ dyn/cm^) 


14.54 


15.74 


18.58 


12.37 


Mg/Mq 


1.4 


1.4 


1.4 


1.4 


{Ma - Mg)/Mq 


0.173 


0.162 


0.169 


0.203 


R (km) 


11.37 


11.15 


10.86 


11.56 


Ac (km) 


1.08 


1.06 


0.88 


1.12 


1(10^4 g ^j^2^ 


13.34 


12.80 


12.30 


13.97 


z 


0.254 


0.261 


0.271 


0.248 
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TABLE VII 



Quantity 


TF90n 


TF90 


TF94 


TF96 


HWW2 


wuu 


ee (MeV/fm^) 


410.6 


430.1 


490.3 


508.1 


451.3 


531.2 


AM/M 


17 


13 


13 


13 


14 


11 


i?eq (km) 


14.50 


13.59 


12.99 


12.82 


13.03 


12.02 


Rp (km) 


10.64 


10.35 


10.21 


10.16 


10.18 


9.82 


log //(g cm^) 


45.0829 


45.0704 


45.0568 


45.0525 


45.0596 


45.0345 


t 


0.082 


0.076 


0.068 


0.066 


0.070 


0.060 




0.610 


0.600 


0.595 


0.593 


0.594 


0.579 




0.2792 


0.2905 


0.2966 


0.2986 


0.2979 


0.3141 


e 


0.68 


0.65 


0.62 


0.61 


0.62 


0.58 


n (km^) 


14.7 


12.6 


10.4 


9.8 


11.0 


8.1 



TABLE VIII 



Quantity 


TF90n 


TF90 


TF94 


TF96 


HWW2 


WUU 


ec (MeV/fm^) 


454.0 


465.8 


531.3 


549.9 


480.9 


560.7 


AM/M 


6 


5 


5 


5 


6 


4 


Rcq (km) 


13.05 


12.47 


12.00 


11.88 


12.07 


11.30 


Rp (km) 


11.74 


11.35 


11.03 


10.94 


11.06 


10.51 


log 1/ (g cm^) 


45.1419 


45.1275 


45.1055 


45.0990 


45.1155 


45.0803 


t 


0.031 


0.029 


0.026 


0.025 


0.027 


0.023 




0.650 


0.635 


0.625 


0.622 


0.625 


0.606 


Zp 


0.2425 


0.2549 


0.2649 


0.2679 


0.2649 


0.2841 


e 


0.44 


0.42 


0.39 


0.39 


0.40 


0.37 


n (km^) 


5.5 


4.8 


3.9 


3.7 


4.2 


3.1 
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TABLE IX 



Pulsar 




Model atmosphere 


T^[K] 


0833-45 


1.1 X 10* 


blackbody 


1.3 X 10^ 


(Vela) 




magnetic H-atmosphere 


7.0tli X 10^ 


0656+14 


1.1 X 10^ 


blackbody 
H-atmosphere 


7.8m X 10^ 

5.3tli X 10^ 


0630+18 


3.2 X 10^ 


blackbody 


5.2 + 3.0 X 10^ 


(Geminga) 




H-atmosphere 


1.7+1.0 X 10^ 


1055-52 


5.4 X 10^ 


blackbody 


7.91^:^ X 10^ 
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Figure captions 



Fig. 1. EOSs of symmetric infinite nuclear matter. Compared are: TF, 
WUU and WUT, BHF2, BHF3, RB, and HWWl. (The abbreviations 
are explained in Table II.) 

Fig. 2. Comparison of pure neutron matter EOSs. 

Fig. 3. Enlargement of the EOSs of neutron matter at low densities. La- 
bels as in Fig. 1. The triangles (HWWl) show the results of DBHF- 
calculations (potential A) performed in the full Dirac space (for details, 
see Ref. [|T], g], H). 

Fig. 4. Comparison of the symmetry energy in the relativistic and nonrela- 
tivistic approaches. 

Fig. 5. Variation of the EOS with composition in the TF-model. Compared 
is pressure versus density of neutron star matter, composed of p, n, e~ 
and n~, with the pure neutron matter case. The influence of muons on 
the EOS is exhibit too. 

Fig. 6. Comparison of EOSs of /3-equilibrated (i.e., n,p, e~ and /i~) neutron 
star matter. 

Fig. 7. High-density behavior of the EOSs of /3-stable neutron star matter. 

Fig. 8. Pressure-density relation of neutron star matter for different models. 

Fig. 9. Comparison of the proton fractions associated with the various dif- 
ferent EOSs studied in this paper. 

Fig. 10. Gravitational mass of static (nonrotating) neutron star sequences 
versus central energy density for different neutron star matter EOSs. 

Fig. 11. Neutron star radius versus mass for TF and the relativistic HWW2- 
EOS. Shown are sequences of stars rotating at their general relativistic 
Kepler frequencies, and at zero frequency. 

Fig. 12. Gravitational neutron star mass versus central star energy density 
for TF and the relativistic HWW2-E0S. Shown are sequences of stars 
rotating at their Kepler frequencies, and at zero frequency. 
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Fig. 13. Mass density profiles of a 1.4 Mq neutron star computed for several 
different EOSs. 



Fig. 14. Limiting rotational periods (P = 27r/i7) of neutron star sequences 
versus star mass. Shown are the Kepler periods and the limiting periods 
set by the emission of gravity waves from young neutron stars (lO^^X) 
as well as old ones (lO^-fC). The underlying EOSs are TP and the 
relativistic HWW2. The shaded area covers the range of observed 
periods and masses (see text). 

Fig. 15. Cooling of neutron stars constructed for TP, WUU and HWW2, 
with and without inclusion of superfiuidity. The surface tempera- 
tures obtained with a blackbody- (magnetic hydrogen) atmosphere are 
marked by solid (short-dashed) error bars (see Table IX). The uncer- 
tainty in the pulsar's age is indicated by the error bar in the lower right 



corner 39 
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FIGURE 1 
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FIGURE 2 
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FIGURE 3 




FIGURE 4 



200.0 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

p (1/fm') 



25 



FIGURE 5 
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FIGURE 8 
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FIGURE 9 
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FIGURE 10 
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FIGURE 11 
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FIGURE 13 
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